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Abstract 

The long-standing problem of constructing a potential from mixed scattering 
data is discussed. We first consider the fixed-£ inverse scattering problem. We show 
that the zeros of the regular solution of the Schrodinger equation, rn{E) which are 
monotonic functions of the energy, determine a unique potential when the domain 
of energy is such that the r„(£')'s range from zero to infinity. The latter method 
is applied to the domain {E > Eq,£ = Iq} U {E = EQ,i > Iq} for which the zeros 
of the regular solution are monotonic in both parts of the domain and still range 
from zero to infinity. Our analysis suggests that a unique potential can be obtained 
from the mixed scattering data {5{io, k), k > /cq} U {5{i, ko),£ > Iq} provided that 
certain integrability conditions required for the fixed ^-problem, are fulfilled. The 
uniqueness is demonstrated using the JWKB approximation. 



* e-mail: lassaut@ipno.in2p3.fr. 
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1 Introduction 



Three-dimensional (3-D) inverse quantum scattering problem amounts to reconstruct 
some interaction potential from a set of experimentally accessible scattering data. For 
spherically symmetric potential^] V{r), the case that will be considered throughout this 
paper, it is known [1] that the 3-D problem (— V{r))ilj{E,r) = Etp{E,r) {Af is the 
3-D Laplacian, E is the energy) can be reduced, thanks to the partial wave decomposition 



tjj{E,r) = X]/gn Sm=-f ^'^^f'^^ ^e,m{d, 0) (the Yi^^s are the spherical harmonic functions) 
to the study of the radial Schrodinger equation with centrifugal barrier term depending 
on the angular momentum i. This latteio is given by 

dr^ r 



i:-,+E- V{r) - S^:^^^ tlZ ) ^^^E, r) = (1.1) 



where ip£{E,r) is called the regular solution which is uniquely defined, as usual [I112], by 
the Cauchy condition \imr^o'^e{E,r)r~^~^ = 1. It behaves for positive values of E as 
^pi oc sin(fcr — £tt/2 + 6{£, k)) when r — > oo (fc = \fE, see second footnote), provided that 
V{r) satisfies the integrability condition [2j 

+ 00 fOO 

\V{r)\dr <oo, b > 0, / r\V{r)\dr < oo (1.2) 



Here, the fc)'s are the phase shifts, the relevant scattering data for the present paper. 
Notice that the whole scheme can be extended to real values for i verifying 2i + l > [Ij. 

Basically, approaches to 3-D inverse scattering can be classified in two categories [U, E] , 
namely, those using scattering data for fixed values of either the energy E or the angular 
momentum i. In the first category, the so-called fixed-i? inverse scattering problem, 
Loeffel [3j obtained theorems predicting a unique potential from the knowledge of the 
phase-shifts 6{i, k) at a specific energy E = k"^, for all (non-discrete) non negative values 
ofA = £+ l/2. When the set of data is reduced to discrete values of A = £ + 1/2 for 
non-negative integer i, Carlson's theorem |1] predicts a unique potential V{r), provided 
that this latter belongs to a suitable class [21 E]- The Newton series permits one to 
construct the potential V{r) from the knowledge of the discrete partial wave scattering 
data [3J. In the second category, the so-called fixed-£ inverse scattering problem, the 
potential V{r) that satisfies (11. 2p can be constructed from the phase-shifts 6{i, k) given 
for all momenta k G (0, oo) and from the discrete spectrum data (bound state energies 
and the corresponding normalization constants). For details, see pj- 

Apart from the above two distinct approaches there exists also a possibility of con- 
structing a potential starting from input scattering data that are partly i^^-dependent and 
partly ^-dependent. The idea of exploiting simultaneously E- and ^-dependent scattering 
data was first explored by Grosse and Martin [5] for the construction of confining po- 
tentials. They conjectured that the knowledge of the ground state energies for all 
non-negative integers i, allows one to recover the potential in a unique way. The prob- 
lem of reconstructing potentials from bound states has been studied numerically in Ref. 



^Notations: f = (r sin cos (p, r sin sin 0, r cos 0),?'>O, 0<6'<7r, O<0< 2tt. 

^The potential v{r) is normalised such that ^V{r) — v{r) where m is the mass for the particle. We 
use throughout this paper the units [2\ h = 2m = 1; then, V{r) ~ v{r) and the energy E verifies E = k"^ 
where k is the wave number. 
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[6]. In [7j, Rudyak and Zakhariev have constructed potentials from E- and ^-dependent 
data, however assuming that aE + bi{i + 1) is a constant. In this case, the corresponding 
analysis is based on an extension of the Newton's method. 

In this paper, we construct a potential obtained from input scattering data that are 
partly E- and ^-dependent, when no extension of the Newton's method is available. Let 
us summarize the results of our study. First we show that, in the fixed-£ problem, the 
knowledge of the zeros of the regular solution the Schrodinger equation permits one to 
determine the potential. More precisely, we show that the zeros of the regular solution, 
denoted hereafter by rn{E), n > 1, which are known to be monotonic functions of the 
energy, whatever the behavior of the potential may be [8], determine a unique potential 
provided the domain of the energy is such that rn{E) range from to oo. Note that the 
reconstruction of the potential from the zeros of the regular solution enters the category 
of inverse nodal problems introduced by Hald and McLaughlin [9J. As a second result, 
we show that piecewise constant potentials can be constructed from only one single line 
of zeros still ranging from to oo, say rno{E), for some fixed uq. As a third result of 
this paper, we prove uniqueness theorems, stating basically that two different potentials 
satisfying f lL2p cannot have a common line of zeros ranging from to oo. In the second 
part of the analysis, we turn to the mixed £ and E problem and apply the above results to 
the domain {E > Eq, i = £o}U{£' = Eq, £ > Eq} for which the zeros of the regular solution 
for the Schrodinger equation are monotonic functions on both parts of this domain and 
still range from zero to infinity. This domain cannot be dealt with (any extension of) 
the Newton's method. The first part of the analysis suggests that a unique i- and E- 
independent potential could be obtained from the set of mixed scattering data given by 
{6{io,k),k > ko} U {6{i,ko),i > io}. As the last result of this paper, we show that 
this is indeed true within the semi-classical JWKB approximation, provided there is only 
one turning point. Note that the uniqueness theorem 1 still applies for the mixed data, 
leading to the construction of piecewise constant potentials from the zeros of the regular 
solution of the Schrodinger equation. General details on the JWKB approximation and 
its full applicability can be found e.g in [IQl while specific applications to the Schrodinger 
equation may be found e.g in IH [2] . Basically, it provides a way to determine a semi- 
classical expansion (in h) of the solution of the Schrodinger equation and is valid whenever 
the potential changes slowly during an oscillation of this solution. 

Our analysis has natural applications in heavy ion physics as well as general properties 
of scattering processes although applications in other areas, such as acoustics and/or 
geophysics seem likely to be possible. In heavy ion (elastic) scattering physics, a natural 
application of the present paper concerns the problem of discrete ambiguities encountered 
in optical model analysis [H] and would at least provide a new insight on past works of 
this problem p^[T3]. We have undertaken the corresponding study which will be reported 
elsewhere [13]. The present analysis can be easily extended to more complicated set of 
(mixed) scattering data. In particular, we emphasize that the most noticeable feature 
of our approach is that it can be applied to inverse scattering problems not reducible to 
the i = 0-fixed inverse scattering problem through a Liouville transformation (defined for 
instance in [2]), as it is the case for the fixed- problem for potentials of finite range [3] 
or for the case treated by Rudyak and Zakhariev [7] also for potentials of finite range. 
Note that, when the inverse scattering problem is reducible to the usual s-wave inverse 
scattering problem, there exists an extension of the Newton's method available for the 
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scattering problem considered. 

The paper is organised as follows. The section [2] involves the whole analysis. In 
subsection 12.11 , we consider the fixed-i problem. In subsection 12.21 we discuss in detail 
the mixed scattering problem. The subsection 12.31 deals with the JWKB approximation 
applied to the mixed scattering data. In section [3], we discuss the results and conclude. 

2 Formalism 

2.1 The fixed-£ inverse scheme 

We assume that the potential fulfills the integrability condition fll.2p . Let us consider 
the regular solution of the Schrodinger equation ip^{E,r) for fixed positive energy E and 
fixed a function of r > 0. It can be realized that ipe{E,r) and ilj'^{E,r) cannot vanish 
at the same value(s) of r for 2£ + 1 > (except for r = 0) and ^ > as shown in [15] 
so that ipi{E,r) can only have simple zeros (except at r = 0). Owing to the asymptotic 
properties of ipAE^r) and ^(E,r) as r — > oo, it can be further realized that the set of 

zeros is countable. According to these observations, these zeros (the origin r = being 
excluded) can then be naturally ordered starting from the smallest value that will be 
denoted by ri{£,E). A generic zero of the regular solution will be denoted by rn{i,E), 
n e N, n > 1. This defines the numbering of the zeros that is used in the following analysis. 
Notice that, thanks to the fact that all the zeros are simple, one has the following (strict) 
ordering < ri(£, E) < r2{i, E) < ■ ■ ■ < rn{i, E) < ■ ■ ■ . 

For any potential, the zeros satisfy the monotonicity properties, i) E = k"^ rn{i, E) 
is a decreasing function, as has been shown by Sturm in 1830 's [H] and ii) (i ^ rn{i, E) is 
an increasing function. This can be easily shown. More precisely, we have 



Note that, as mentioned above, ipe{E, r) and {dip£/dr){E, r) cannot vanish simultaneously 
for 2i + 1 > 0, except for r = and i > [T3], so that the denominators in fl2.ip and 
(12. 2p cannot vanish. 

For potentials satisfying ( fll.2p ) the function E h-* rn{i,E) is such that rn{i,E) 
for oo [16j. This is a consequence of equation (1.5.6) of [2]. In the absence of bound 
states rn{i, -E) — > oo for ^ 0. 

For negative or zero values of E, the regular solution of the Schrodinger equation has 
zeros provided that the potential has bound states, and that E > Ei, where Ei < 
denotes the ground state energj{§. The number of zeros is finite for values of the energy 
E such that Ei < E < 0. For details see [16]. When the potential has N bound states. 
El < E2... < E]\f, if n < N, we have r„(£, E) ^ 00 for E ^ En [IS], whereas for n > N, 
rn{i, E) 00 a.s E ^ 0. This is illustrated in the Fig. [T] where we have drawn the first 
four zeros of the regular solution for a Bargmann transparent potential of Ref . [TT] , which 

■^We adopt this notation to denote the ground state by Ei instead of the traditionally used Eq in order 
to be consistent with the meaning given to n, namely to denote zeros of the wave function 
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Figure 1: Zeros rn{E), n = 1,2,3, and 4, of the regular s-wave solution for a Bargmann 
potential (for more details see text). 

has a bound state at the energy E = —1 in units. All zeros, but the first one, go 
to infinity as E goes to 0. The first zero, however, goes to infinity when E ^ Ei = —1. 
In the fixed-^ inversion, if all the zeros are assumed to be known, i.e. E \—>- rn{i,E) is 
known and r„(£, E) describes the entire interval [0, oo[, then the potential Ve{r) is uniquely 
determined. This can be easily checked. For this we consider the Sturm-Liouville problem 
on [0,i?], i.e., the equation 

^eir) +[e- Ve{R - r) - Mr) = 0, (2.3) 

coupled with the Dirichlet conditions 

^^^(0) = ^,{R) = . 

The spectral data are the eigenvalues -E*'s, such that R = rn{i, E*), and the normalization 
constants \IS\ 

The latter are given by the positive values 

Pn = -^rn{i,E)\E=E^. 

The potential, assumed to be square integrable on [0, R] is uniquely determined on [0, R] 
by the spectral data {E*, p„} [19\. The technique of constructing V from this set of data 
is well-known |2l HI [20l [2l] . 

If we consider now the Sturm Liouville problem on [0, R] consisting of the equation 

^;'(r) +(^E- Ve{r) - Mr) = 0, (2.4) 

coupled with the Dirichlet conditions 



MO) = MR) = 



we recover the Gel'fand-Levitan approach at the hmit i? — > cxo [2], [22] . 

For the purpose of using mixed, E- and ^-dependent data, we start with the following 
remark. For a fixed i, the nth zero of the regular solution can be considered as a function 
of energy r„(£, E), which is monotonic and such that ipe{E, rn{i, E)) = 0. It defines a line 
of zeros as we mentioned in the introduction. Moreover E ^-^ rn{i,E) admits an inverse 
function r ^-^ En/{r) which is also monotonic and is the inverse of this line of zeros. For 
example, consider the potential V = in the S-wave {i = 0). The regular solution is 
proportional to sm(\/Er). The lines of zeros are given by r„(0,-E') = mr/^/E and the 
inverses of these lines are E'^.o = n^vr^/r^. 

Suppose that the inverses r h-* En.e, of all lines of zeros n > 1, are known for r running 
from to oo. This implies that E* = En/{R) is also known and that the p„ are given by 

— = --J-En,eir)\r=R 

Pn dr 

We can then construct the desired potential, assumed to be locally square integrable on 
[0, R]. If only the inverse r i— > E^/ of a single line of zeros is known, no method is available 
to recover the potential, except in the special case of piecewise constant potentials. What 
we can show is a uniqueness property, if En^^i{r) is known for all positive r (see the 
theorems given below). 

In the special case of piecewise constant potentials, having discontinuities at values of 
r = Oj, j = 1, . . . ,jmax, and being zero for r > aj^^^, it is easy to show that there is a one 
to one correspondence between the discontinuities of the third derivative of i?„j/(r) with 
respect to r and the discontinuities of V. This suggests the following lemmas: 
Lemma 1: For a piecewise constant potential, the knowledge of a single line of zeros 
r{E) = rn{i,E) (n fixed) allows the reconstruction of the potential in a unique way pro- 
vided that the latter function has values ranging from zero to infinity when E e]0, +oo[. 

The following variant of the lemma 1, although useless for the mixed set of data con- 
sidered in the next subsection, works also for potentials having bound states. Denoting 
r \^ E{r) the inverse function of the function E i-^ r{E) which is the single line of zeros 
r{E) considered, we state that 

Lemma 2: For a piecewise constant potential, the knowledge of the inverse of a single 
line of zeros E{r) = En/{r) allows the reconstruction of the potential in a unique way 
provided that the latter function of the variable r is known on the whole interval ]0, -|-cxd[. 
The proof uses the arguments given above, in particular the monotonicity property of the 
line with respect to the energy, and the third derivative of 

ip{E{r),r) = 0. 

To be more specific, if E'" has a discontinuity at r = a then 

^(«*) - -^i"-) = -^^i") [f - n--)] ■ (2-5) 

This is equivalent to the relation 

{Et)-^^{E~) = 2l^—{E^)j [V{r{E:))-V{r{E~))] (2.6) 



dE^ 
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Since the potential is zero for r > aj^^^, the relation (12.51) allows us to reconstruct the 
potential between aj^^^ and The procedure can be repeated at each aj, and the 

potential is obtained at r 7^ aj by summing the successive values at each discontinuity 
appearing beyond r. 

As an illustration consider the potential 

{-2 r < 2 
-1 2<r<3 (2.7) 
r > 3 

In Fig. [21 we have drawn the function r i— > —E"'{r)/{2E'{r)), related to the inverse E{r) 
of the first line of zeros for the S-wave regular solution of the Schrodinger equation with 
the potential (12. 7p . Clearly, the discontinuities of r —E"'{r)/{2E'{r)) happen at the 




1 1.5 2 2.5 3 3.5 4 

r 



Figure 2: The structure of the function —E"'{r)/{2E'{r)) related to the inverse E{r) of 
the first line of zeros for the S-wave regular solution corresponding to the potential (12. 7p . 

points where V has discontinuities, namely r = 2 and r = 3, and the equation (12.51) is 
satisfied. We know that the potential is zero beyond r = 3. So we have = ^(3"^). 
From the curve the discontinuity ^(3"^) — ^(3~) is equal to 1 then ^(3~) = —1. As V is 
piecewise constant, it is constant on ]2,3[ and we have V{2~^) = V^(3^) = —1. From the 
curve the discontinuity V{2~^) — V{2~) is equal to 1 so that V{2^) = —2. We thus recover 
the potential (12. 7p . 

This method cannot be applied to a potential defined by a continuous function. Nev- 
ertheless, for such potentials, the following uniqueness theorems hold: 
Theorem 1: Consider two potentials V\ and V2 satisfying i\1.2\] and locally constant in 
the vicinity of zero. For the fixed-i problem, assume that both potentials have no bound 
states. Suppose that two integers ni and n2 exist such that the ni-th line of zeros for the 
regular solution for V\ coincides with the n2-th line of zeros for the regular solution for V2, 
(rnj^{£, E) = Tn^i^, E) ) with both lines having values ranging from to 00 as the energy E 
varies in ]0, oo[. Then V\ = V2. 

Its variant, still useless for our mixed set of data, but working for potentials having 
bound state yields: 
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Theorem 2: Consider two potentials V\ and V2 satisfying (11.21) and locally constant in 
the vicinity of zero. For the fixed-i problem, assume that two integers ui and n2 exist such 
that the inverse of the ni-th line of zeros for the regular solution for Vi coincides with 
the inverse of the n2-th line of zeros for the regular solution for V2 on the whole interval 
]0,+oo[, ((Vr > 0)E„,,,(r) = ^„„^(r)). Then Vi = V2. 

Let us prove the theorem 1. The proof of its variant, the theorem 2, quite similar to that 
of theorem 1, is not reported here. 

First let ilJi{E,r) and ip2{E,r) (we suppress the index i for convenience) be the regular 
solutions of the Schrodinger equation for potentials Vi and V2, respectively, constrained 
by the Cauchy conditions linv^o "^iiE, r)r~^^^ = 1, i = 1, 2 in the vicinity of zero. For an 
S'-wave, this is equivalent to ipl{E,0) = 1, i = 1,2, the prime being the derivative with 
respect to r. Setting = Vi — V2 and noting that the Wronskian is zero for r = we 
obtain 

W{ij,,iJ2) = ij[{E,r)ij2{E,r)-ij,{E,r)ij',{E,r) 

dr' AV{r')^i{E,r')tlj2{E,r') . (2.8) 



For r = r{E) (we adopt the convention r„j(£, E) = r„2(^, E) = r{E)) we have 

rr{E} 

{WE) / dr' AV{r')iPi{E,r')iP2{E,r') = . (2.9) 



The latter integral exists since ipi, i = 1,2 are continuous and both potential Vi and V2 
satisfy fll.2p ) and are locally constant in the vicinity of zero. 

Differentiating twice Eq. (12.91) with respect to E, which is feasible as the first and the 
second derivative of ipi{E, r)ip2{E, r) with respect to the energy E is a function continuous 
of the variable r [23j, and taking into account that il)i{E,r{E)) = 0, i = 1,2 for every 
value of E', we obtain 

pr{E) r.2 

(WE) dr' AV{r') [^,{E, r') ME, r')] = (2.10) 

We recall that, due to its monotonicity, E 1— > r{E) admits an inverse r 1— »• E{r). We define 
the Kernel 

K{r,r') = ^[MEir),r')MEir),r')] , r' < r . (2.11) 
The diagonal part, K{r,r), equal to 

d , , ^ ^ d , , ^ , „ f dr \ d , , „ , d 



K{r,r) = 2—ME{r),r) —ME{r),r) = 2 J —ME{r),r) —ME{r),r) . 



(2.12) 

is never zero, since ipi and -^"ipi {i = 1, 2) do not vanish simultaneously, except for r = 
and for £ > [15]. 

Introducing the definition (12.111) . the equation fl2.10p is equivalent to 

(Vr > 0) / dr'K{r,r') AV{r') = . (2.13) 

Jo 
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Since the potential Vi (V2) is locally constant in the vicinity of zero, there exist Vi and 
ei ( V2 and €2) such that 

(Vr<ei) V^{r) = V, 

(Vr<e2) V2ir) = % (2.14) 

In the vicinity of zero, the ni-th line of zeros of the regular solution for Vi is given by 

r„,(£,E) = -,fe^ (2.15) 
\/E~Vi 

where ji^m is the ni-th zero, except for the origin, of the function x ^ je{x) involving the 
spherical Bessel function je{x) [21]. 

This implies that 

Taking into account fl2.14p and (12.161) as well as the identity = we see that for 
r < e = inf (ei, 62) we have Vi — V2 = 0. Therefore AV{r) is zero for r < e. 

Differentiating (I2.13P with respect to r, ( possible since the derivative r' \—>- Kr{r,r') 
is continuous) and taking into account that AV is zero for r < e, leads to the Volterra 
equation 



where 



dr'Ki{r,r')AV{r') = (2.17) 

K (r r') BK 
i^i(r,r') = 9pY, ir,(r,r') = ^(r,r') (2.18) 

K (r, r) or 

Let us consider fl2.17p for r G [e, B] . The Kernel Ki is bounded and continuous. This 
is due to the fact that K{r,r) is continuous and never vanishes on the compact interval 
[e,R]. So \K{r,r)\ is bounded from below by a strictly positive value. Thus Eq. (12.171) 
has a unique solution for r G [e, R] [23], which reads AV{r) = 0. By increasing R we can 
conclude that AV{r) = for every r. This terminates the proof. Note that the fact that 
Ki{r,r) behaves like 1/r in the vicinity of zero forbids us to conclude that (I2.17P has a 
unique solution for e = [25]. This is the reason why the potentials have been assumed to 
be locally constant in the vicinity of zero. 



2.2 Inversion scheme from mixed data 

Consider the set {E > Eq, i = £0} U {E = Eq, i > io}. Similarly to Sect. 12.11 we may de- 
fine lines of zeros in which the n-th line of zeros r„(£, E) of the regular solution describes a 
line formed of two parts. In the first part, the zeros r„(£o, E) range from r = 0, {E = 00), 
to Tq = r(£o, Eq), [E = Eq), as the energy E varies from 00 to Eq {Eq being fixed); in the 
second part, r„(£, £"0) has values ranging from tq = r(£o, Eq), {i = Eq), to 00, (£ = 00) as 
£ ranges from Eq to 00 {Eq being fixed). This has been verified for potentials satisfying 
the integrability conditions (ll.2p [2\. The monotonicity property required in the above 
lemma 1 and theorem 1 is preserved in both domains. 
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The lemma 1 is still valid for piecewise constant potentials. In the special case where 
a discontinuity v = V^Tq) — V^Tq) appears at the junction point between both parts of 
the domain, i.e. at tq = r„(£o,-Eo), the value of v can be determined as follows. One 
first adds v to the bare value Vq of the potential at the origin, obtained by adding all 
discontinuities (except v). Then, for r close to the origin, the zeros are given by 

rn{io,E) = ^=M^^= (2.19) 

and the value of v can be obtained from the relation 

dlnrn{io,E) 1 1 



dE 2 E-Vn-v 



(2.20) 



To demonstrate the applicability of the theorem 1, we consider potentials satisfying the 
integrability conditions (11. 2p and locally constant in the vicinity of zero. Using the Wron- 
skian relation (12. 8p . we have to show that, if 

/ dr' A\/(r')^i(4, E, r')^2(^o, E, r') = for E > , (2.21) 



f 

Jo 



and 

'■r{£,Eo) 

dr' AV{r')^i{iEo, r')^2(^, ^o, r') = for £ > 4 , (2.22) 

'0 

then Vi = V2. (Here we have put as previously = r„2 = r.) 

This can be proved in two steps. First, for r < tq = r(£o,-E'o), Eq. (I2.2ip is zero for 
every E > Eq. Consequently, Vi = V2 for r < tq when Vi, V2 are locally constant in the 
vicinity of zero. 

Secondly, taking into account the above argument, Eq. (I2.22p can be rewritten 



rrii,Eo) 

/ dr'AV{r')M^Eo, r')^^, Eq, r') = (2.23) 



Differentiating twice with respect to i, which is feasible as the first and the second deriva- 
tive of ipi{^Eo,r)ip2{^, EQ,r) with respect to £ is continuous with respect to the variable 
r [23], and inverting the monotonic function i r(£, E'o), denoted hj r i{r,Eo), we 
obtain 

(Vr>ro) f dr' AV{r')K{r,r') = (2.24) 

J TO 

(Vr' < r) K{r, r') = ^[V'i(^(r, i?o), i^o, r')V'2(^(r, i?o), i^o, r')] . 
The diagonal part K{r,r) is equal to 

K{r,r) = 2-MKr,Eo),Eo,r)-MKr,Eo),Eo,r) 



2 1:7^) Q^M^{r,Eo),Eo,r)—M^{r,Eo),Eo,r), (2.25) 
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and does not vanish for r > tq > 0. 

Differentiating Eq. fl2.24p with respect to r, (possible as r' i— > Kr{r,r') is continuous) 
we obtain the Volterra equation 

(Vr>ro) AV{r)+ dr'Ki{r,r')AV{r') = (2.26) 

Jro 

where Ki is given by f l2.f8p . Consider now Eq. (12.26 1 ) for r G [^c-R]. Since K{r,r) in Eq. 
(12.251) is continuous on [ro,-R], its absolute value reaches its strictly positive minimum 
m. The Kernel Ki{r,r') is then bounded and continuous and the equation (I2.26P has a 
unique solution AV{r) = on [ro,-R]. Increasing R, we conclude that AV{r) = for 
every r > ro and then Vi = V2 on the whole half axis. 

We have shown that a single line of zeros, which, for the data considered, always has 
values ranging from zero to infinity and moreover is monotonic, determines the potential 
uniquely. The remaining question is to examine whether the set of mixed data 

{6{i = io,k) ke[ko,+oo[}U{6{i,ko)ie[io,+oo[} (2.27) 

associated to the set {i = £o,k E [kQ,oo[} U {k = kQ, i E [£o,oo[}, determines a line of 
zeros, and thus the potential, in a unique way - which is suggested by the analogy with 
the £-fixed problem. In the absence of bound states, all lines of zeros E 1— > rn{i,E), 
monotonic with respect to E, range from zero {E infinite) to infinity {E = 0) when E 
has values ranging from infinity to zero. In this case, we know that the potential, when it 
satisfies (II. 2p . is recovered in a unique way, given the phase-shifts 6{i, k) for all > 0. The 
condition (11.21) excludes all pathologies; for example, potentials behaving asymptotically 
like 1/r^, encountered in particular in the presence of a zero energy bound state, or ghost 
components in the Jost function [TTj. With the mixed data we are in the same situation, 
namely all the lines of zeros are monotonic and range from zero {E infinite, ^ = £0) to 
infinity {E = Eq, i infinite). Thus we expect that (I2.27P is associated with a unique 
potential decreasing faster than 1/r^ at infinity. We cannot prove this in the general case, 
but we can investigate the problem in a JWKB approach. 

Note that a way to treat the general case is to set a self-consistent procedure. First 
we note that, in order to calculate the phase-shift say 6{i,k), there is no need to know 
the potential for distances smaller than r„(£, E), whatever is n. It is sufficient to consider 
the Schrodinger equation on the interval [r„(£, E),oo[ because only the ratio ipe/ip'^, where 
the prime denotes the derivative with respect to r, has to be known to determine the 
phase-shift. In this respect, we first determine the potential in terms of the phase-shift 
6{i, kQ){6{io, k)) and the n-th {n fixed for all the procedure) zero of the regular solution 
by solving the Schrodinger equation on [r„(£, /cq), C)o[ ([r„(^o, fc), C)o[) because the phase- 
shift does not depend on the potential V{r) for r < r„(£, A;o)(r„(£o, fc))- The algebraic 
construction of the potential is guaranteed by the fact that the monotonicity property of 
the n-th line of zeros of the regular solution in terms of the energy does not depend on 
the potential. Then we have to determine the zeros. For instance take the zeros of the 
free solution ^ = then determine V then recalculate the zeros etc.. 

We first consider results based on the Born approximation. In the seventies Reignier 
used a Born approximation of the scattering amplitude to show that the knowledge 
of the phase-shift at a fixed energy, Eq = k^ say, for each integer i is equivalent to the 
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knowledge of the Fourier sine transform of the potential rV{r), 

POO 

9(1) = / sm{qr)rV {r)dr , (2.28) 
Jo 

for q < 2kQ. The scattering amplitude is determined from the phase-shifts at fixed energy 
5{i, ko) for i = 0,1,2,..., E = Eo = kl 

Generally, the integral is assumed to be zero for q > 2ko [21 [12] leading to potentials 

rV{r) = — / sm{qr) g{q) dq , (2.29) 
Jo 

such that rV{r) is an entire function of r of order 1. Other extensions of g{q) are studied 
in [27]. 

More recently Habashy and Wolf [2H] have studied the reconstruction of a potential, 
having compact support and spherical symmetry, from its 3D-Fourier transform through- 
out the Ewald limiting sphere, \ko{u — u')\ where u, u' are normalized to unity and take all 
possible directions. This is equivalent to the Fourier transform of the potential for values 
of k such that \k\ < 2ko. In this case, the inverse sine transform rV{r) coming from Eq. 
(12.291) is given in terms of the potential V{r) by 

ri/(r)= rdr'r(r')^^^^i^^^^^^^^, -R<r<R, (2.30) 
J-R 7r(r - r') 

where R denotes the support of V. For negative values of r, the authors of [2H] have put 
V{—r) = V{r). The latter equation can be solved numerically, thanks to the spectrum 
of the kernel involved in its right hand side, namely, to the z/„ g]0, 1[, ipn G 1]) 
(depending on ko) such that 

, ,sin[2A;o(a; — x')l , , , , x , 

dx' ^ . °^ i)^{x') = M^) , -1 < X < 1 . 2.31 

_1 TC{X - X') 

Coming back to Eq. (12.291) . a possible way to extend our knowledge of g{q) beyond 2kQ is 
to take the Born approximation for the missing phase shifts 6{i = 0, k) for k > ko. This 
is given by 

POO 

/ sin(A;r)2 V{r)dr = -k5{i = 0,k) . (2.32) 
Jo 

The derivative with respect to k yields 

g{q)= r sm{qr)rV{r)dr = -^^^^^^^^^^ , ^q = 2k > 2ko . (2.33) 

This implies that q{q), known for q < 2ko, is now known for every positive q, and that 
V{r) is uniquely given by 

rV{r) = — I sin(gr) (7(g) dg . (2.34) 
TT Jo 

Consequently, the knowledge of {5{C,ko),C- G Af} U {5{C = 0,k),k > ko} allows us to 
recover the potential in Born approximation if ko is sufficiently high. 
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2.3 The JWKB approximation 

In what follows, we examine to what extent the mixed data fl2.27p allow us to recover 
the potential in the semi-classical JWKB approximation. General details on the JWKB 
approximation and its full applicability can be found e.g in [10] while specific applications 
to the Schrodinger equation may be found e.g in [H [2]. It has also been extensively used 
in the fixed energy inverse scattering (see, for example, [291 [I2l [IS [30] and references 
therein). Basically, we recall that this semi-classical approximation provides a way to 
determine a semi-classical expansion (in h) of the solution of the Schrodinger equation 
and is valid whenever the potential changes slowly during an oscillation of this solution. 
In the JWKB approximation, the phase-shift is given by 



6{i, k) 



lim ( / K{\k,r')dr' - [ Kf,ee{\k,r'))dr' (2.35) 

V Jr{\,k) Jri,,,{X,k) J 



(here = fc^ in 1/L^ imits and X = i + 1/2). In the absence of a Coulomb component, 
we have 



K{\, k, r) = ^k^ - V{r) - AVr2 , Ki,^e{\ k, r) = - X^/r^ . (2.36) 

In f l2.35p . r(A, k) is the turning point relative to the function K considered, assumed here 
to be unique for the sake of simplicity. More precisely r(A, k) is a solution of the equation 
K{X,k,r) = 0. 

For the part of the spectrum 6{X, ko), X > Xq we use the results of Sabatier [12\ and 
Cuer [13] for a single turning point. We assume that (A; = ko being fixed) for every A, the 
equation 

kl - Vir) - AVr^ = (2.37) 

has a unique solution denoted by r(A, ko). This happens when the function r i— > g{r) = 
r'^ikP' — V{r)) is monotonic. For V = Q the equation fl2.37p has a unique solution, denoted 
by ?"free(A, k) = X/k. In [121 [13] the phase-shift is given by: 



SUM- I dA'yA^-,n(^-^j (2.38) 

where V{r) is assumed to be differentiable to ensure the derivative of the turning point 
with respect to A exists. The function 

A^^(A) = A^lnf^%^') (2.39) 

dA Vnree(A,A;o)y 

is assumed to be Lebesgue integrable. 
From Eq. f l2.38p we have 

the derivation being possible for A 7^ as the integrand is dominated on ]A, oo[ by the 
integrable function A' ^ |^(A')|/(AVA'2 - A^). 
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Using the Lebesgue-Fubini theorem [31], [32] which works when dS{i, ko)/dX G L^ijV) 
and bounded, as made for the Abel transform, the turning point is given in terms of the 
phase-shifts by 

Vrfree(A,fco)/ dA' VA'2 - A2 

provided that the logarithm is zero for A infinite. The above hypothesis that d(5(£, fco)/dA 
is bounded avoids pathologies in the integration in the vicinity of A. For a single turning 
point we recover the result of Loeffel: the knowledge of 6{i,k) for every X = i + 1/2 
positive allows us to recover the potential. Indeed, for A infinite the logarithm tends to 
zero and r(A, ko) tends to infinity. When A tends to zero the equation (12.371) has as unique 
solution r = when lim^-^o V{r)r'^ = 0. 

For r > vq = r(Ao, ko), the potential is given by 



rfree(A, kg) 

r(A, fco) 



(2.42) 



Note that the behavior of the potential for r(A, ko) or A infinite is related to the behavior 
of the phase-shift S{i, ko) for A infinite. 

For r < ro, we introduce the turning point r = r(Ao, k) as the solution of 

e-Vir)-^ = 0. (2.43) 

We assume that the latter equation has a unique solution. This happens when r ^-^ 
h{r) = V{r) + Ag/r^ is monotonic, for example for > or for V attractive but "weak" 
enough when compared to the centrifugal barrier, in the sense that ||r^l^'(r)||oo < Aq . 
The aforementioned function h is positive and infinite at the origin and zero for r = oo. 
Thus, we are sure to have at least one solution. The monotonicity of h implies the 
uniqueness of the solution of the equation h{r) = k^, denoted by r(Ao, k). It implies also 
the monotonicity of t— > r{XQ,k), which is monotonically decreasing with respect to k. 
For k infinite r(Ao, k) 0. 

For V" = Eq. f l2.43p has a unique solution, denoted by rfree(Ao, /c) = Xo/k. The 
equation (12.351) is rewritten as 



Siio,k) 



Vk^ - A (^riXo, k') - rfrec(Ao, A;')) dk' . (2.44) 

When V is differentiable, the derivative of the turning point with respect to k exists. We 
assume that the latter is locally integrable. The above equation is derived with respect 
to k and thus 

'5U4,fc) = ^m^.k) = I' A (^riXo, k') - rfree(Ao, k')^ dk' . (2.45) 

Using the procedure of the inverse Abel transform [231 E3| we obtain 



u 



^{u^ -k^){k'^ ~u^) 

(2.46) 
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using the Lebesgue-Fubini theorem [3T],[32] which works when (5^(£o5 k) is Lebesgue inte- 
grable and bounded. Once the integration over u is performed we obtain 

r(Ao, k) - rfree(Ao, k) = - ,f °' , / dfc^ , (2.47) 



vrJo - k'^ 

provided that r(Ao, k) — rfree(Ao, k) is zero for = 0, which happens for finite range 
potential In this latter case, for k high enough, r(Ao, k) = Xo/k and tends to zero for k 
infinite. 

The potential 

V{r{Xo,k)) = e~-^^ 
r(Ao, ky 

can be determined from fl2.47p for < r < R = r(Ao, k^), i.e. for k > k^. However, Eq. 
02.471) requires knowledge of the phase-shift 5(Ao, k) for every value of k, whereas they are 
known only for k > ko- 

Nevertheless, these phase-shifts for k < ko can be determined from the potential 
beyond the distance R, R = r(Ao, ko). Indeed, for k < ko, f l2.35p can be written as 



S{io,k)= r Jk^-v{r{XM))- al^/'nY"^ r r'-^F^^- 

J\(k) V Aofco/fc V ^ ^0 

(2.48) 

Here, A(fc) corresponds to the value of A for which the first square root in the previous 
equation vanishes, i.e it is a solution of 

- kl + 4vT^ = • (2-49) 



r{X,koY 



This implies, for k < ko, that 



with r(A, ko) given by (12.411) in terms of the phase-shifts 5(A, ko), X > Aq. 

3 Conclusion 

In the present work, we have been concerned with a non-standard inverse scattering 
problem, namely with the construction of the potential from scattering data information 
which involves the use of phase-shifts 

{6{io,k),k>ko}U{6{i, ko),i>io} 

corresponding to the domain 

{{io,k),k>ko}U{{e, ko),i>io} 

and without requiring any extension of Newton's method for inverse scattering problems. 
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First, we have shown that, in the fixed-£ problem, the knowledge of the zeros of the 
regular solution the Schrodinger equation permits one to determine the potential: in fact, 
the zeros of the regular solution r„(i?), n > 1, determine a unique potential provided the 
domain of the energy is such that rn{E) has values ranging from to oo. As a second 
result, we have shown that piecewise constant potentials can be constructed from only 
one single line of zeros (say rno{E), for some fixed no) with values still ranging from to 
oo. We also proved uniqueness theorems, stating basically that two different potentials 
satisfying (11.21) cannot have a common line of zeros ranging from to cx). Furthermore, 
we have also considered the mixed i and E problem and applied the above results to the 
domain {E > Eq, i = io} U {E = Eq, £ > Eq} for which the zeros of the regular solution for 
the Schrodinger equation are monotonic functions on both parts of this domain and have 
values still ranging from zero to infinity. This domain cannot be dealt with any extension 
of the Newton's method. A unique i- and i?-independent potential could be obtained 
from the set of mixed scattering data given by {6{io, k), k > ko} U kQ),£ > Iq}. As 
the last result, we have shown that this is indeed true within the semi-classical JWKB 
approximation, provided there is only one turning point. 

Furthermore, we have shown that in Born approximation the following mixed scattering 
data 

{S{i = o,k),ke [ko, +oo[} u {§{£, ko)Je ^^} 

lead to an unique potential, still assumed to be i- and i?- independent, which is the inverse 
Fourier sine transform of a function deduced from the data. 

To conclude, our method, which does not use an extension of Newton's method, takes 
advantage of a set of mixed scattering data. 

A natural application concerns the discrete ambiguities (i.e. the set of different families 
of potentials with practically the same fixed energy phase-shift) encountered in heavy ions 
elastic scattering optical model analyses [11] . These latter ambiguities have been studied 
within the JWKB approximation in Refs. [Ill fT3] . 

Finally we point out that the flexibility of our method, based upon the properties 
of the zeros of the regular solution, provide a way to study inverse scattering problems 
which are not reducible to the i = 0-fixed inverse scattering problem though a Liouville 
transformation. This latter class of problems involves the inverse scattering problem from 
our mixed set of data, for which there is no expected available extension of the Newton's 
method. 
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